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The charge radii and the decay rates of the pion and kaons are calculated, using the relativistic 
equation of motion with a linear potential. Those physical quantities are quite well explained with 
the current quark masses in the case of the pion. It is found that the Van Royen-Weisskopf paradox 
can be cleared out only in the linear potential model by considering the color degree of freedom of 
the quark in the meson. The physical quantities for the kaon are not so satisfactory that it should 
be required to reconsider the Cabibbo-Kobayashi-Maskawa matrix element in the kaon decay. 
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^jQ, Recently the author has attempted to derive a quark confinement potential from quantum chromodynamics (QCD) 
^ 101 . Unfortunately the attempt does not yet succeed, because there is a serious problem to define the Fourier transform 
of a confinement potential which is non-local. However it is natural to think that there should exist a certain confine- 
ment potential in nature, considering hadronic phenomena. So there are some attempts to define it by introducing a 
small parameter HJ^. This letter has the purpose to explain some physical quantities using the potential in order to 
' ' confirm the existence of such a confinement potential in nature. The electromagnetic charge radii and the decay rates 
I , of the pseudoscalar meson are very important physical quantities for understanding the structure of the hadron and 
^ QCD of its theory. Those physical quantities have been calculated many times since the advent of the quark model 
CO However there are some difficulties to describe hadronic phenomena, such as the constituent quark masses 

1^ and the Van Royen-Weisskopf paradox in those literatures. This letter shows that the proper relativistic description 
, explains such difficulties quite reasonably with the current quark masses. The consideration of the color degree of 

• freedom clears out the paradox excellently, which is the difference of the spatial wave function |^'(0)p between the 
^\ , calculations of the charge radius and the decay rate of the pion. In the case of the kaon, however, there is a difhculty to 
0^ ' explain its charge radii and decay rate consistently. The reason is conjectured that the Cabibbo-Kobayashi-Maskawa 

, (CKM) matrix element |KisP is too large so that the kaon decay constant is too small. Therefore it is necessary to 

• scrutinize the CKM matrix element for the kaon decay. 
In order to describe the relativistic feature of the quark in the bound system of the meson, the Klein Gordon 

equation is used to do it rather than the Dirac equation. Since the quark is, of course, a fermion, the Dirac equation 
] should be used to do it. However there are good reasons to adopt the Klein Gordon equation, because the degeneracy 
of spin singlet and triplet states is broken larger than the masses of pseudoscalar mesons so that the spin structure 
of the equation of motion may be less important. Moreover the relative motion of the quark in the meson can be 
assumed to be governed by the effective potential of the spin singlet state. While the electron in the hydrogen atom 



^ J has two spin states to occupy in the ground state regardless of the spin of the proton, the quark has only one spin 
state in the ground state of the meson due to the influence of the other quark, because they have similar masses. This 
is another reason to adopt the Klein Gordon equation, because it is not necessary to describe the motion of the quark 
with a spinor. 

The Klein Gordon equation for a free quark is given by 

(□ + m2)^' = 0. (1) 

Since two body system can be separated into the relative and the center of mass parts, if the equation is modified as 
follows 

-— V?*-(^-mi)* = 

TOl TOl 
1 

Vn^-f— -TO2)* = 0. (2) 

TO2 TO2 
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Using the identity 



—V? + — = -Vj + ^ V% (3) 



where fi is the reduce mass of mi and m2, the sum of the two equations is rewritten as 



-Iv^vl; ^ vIj* - + = 0. (4) 

fi mi + m2 mi TO2 

At the center of mass frame, Vj^'i' — and Pi — P2 = P^, the equation is written simply as 

1 1)2 

--V^*-— * = 0, (5) 
M A* 

where p'^ = Ef — ml = i?! — m| = E"^ — fi^ — p^. The relative energy-momentum and mass are Er = m2E[/{mi + 
m2),Pr — tTL2p'i/{mi + m2), and 11 = m2mi/(mi + m2), where E[ and p'l are the energy and momentum in the rest 
frame of m2. Thus the relative part of the two body system is reduced to the equation of motion of a single particle: 

(a + fi^)^ = 0, (6) 

which satisfies the relative energy-momentum relation. Since the two particles interact with each other with exchanging 
their gauge field, the derivative should be replaced to the covariant derivative as follows 

{D^Df" + fi^)^ ^0. (7) 

where = 5^ — igAf^. A few algebraic effort can make the equation be decomposed as 



w ^ ^ * + ir^^^aAoEr ~ 2gA ■ p, - + g' A')^ = {Er - m)*, (8) 

Er + fi Er + fi, 

where the color index is hidden in the case of QCD. For the test of the validity of the equation, let's approximate it 
for the case of QED i5 « /x, namely, Er = fL + Kr and Kr ^ fi- Thus the equation is expanded as 

- ^V^^- + eAo* - j^A ■ pr-^ + ^AH - ^[{eA^f - {KreAo) + = Kr^^ (9) 

where higher order terms are ignored. If the potential ^0 — ~e/r is inserted, this equation is nothing but the 
Schrodiger equation with the relativistic correction terms in the square bracket which agree to that in the standard 
texts by inserting the relation Kr = —V'^j^fi + eAg in the brace bracket. In the nonrelativistic case of the motion of 
the particle, the potential which make a bound system is the coulomb one. However the potential should be modified 
as the effective potential as insisted in Ref |^ due to the effect of the vacuum polarization of the quark, when the 
motion of the particle reaches to the ultra relativistic region. Hence Eq. (^) can be written as 

-V^^ + 2Er{V^ff.-E^+fi'^)^^0. (10) 

Notice the coefficient 2Er in front of the effective potential which is 2/i in the nonrelativistic case. It seems to 
be the origin to introduce the constituent quark masses in the calculations of low energy physical quantities with 
the Schrodinger equation and other relativistic approaches. For the confinement potential Vcs. = kr, where k = 
2as/37rS/mj under the condition 4my < q^i^^^^ in Ref. [Q, the solution of the differential equation is just the Airy 
function which is 

^{r) = Cfir)e~^^''\ (11) 



where C — ^Jia? j^-K for /(r) = 1 which gives the finite value |5'(0)|2 in the calculations of decay rates. The parameter 
a is calculated as yJ%Erk/9 from the differential equation for the ground state of the quark. Since the coefficient of 
the wave function depends on the energy of the particle instead of its mass, it is natural that the minima for the s- and 
d-wave components of the deuteron function should be shifted toward large momenta in relativistic case as pointed 
out in Ref. Q. In this reference, the potential should be also transformed or rewritten by a scale transformation in 
the momenta, but such a strange thing does not occur in this approach, because the coefficient 2Er plays the role 



2 



sufRciently. Moreover no Lorentz contraction leads to no transformation, because the variable r of the potential is 
perpendicular to the motion of the particle. 

The form factor and the charge radius are defined and calculated for the wavefunction as 



< r2 > = 47r / p(rydr = r(^)(2a)-4/3, (12) 



where Q means the electromagnetic charge of the meson and T designates the gamma function. Since the relative 
distance is not the real charge radius of the pion, the density should be replaced to p[r) — Qu\^u{ru)\'^ + Qd\^ d{rd)\^ 
with = m,dr/ {rriu + md) and = m^r j (rn^ + m^). The charge radii for the pion and kaons are calculated as 

<-^2 ^ 1 .10. 2mg + m^ 9 ,^3 
3^3 ' (mu^mdY^i^E^k' ' 

2 _ 1 10 2ml ( 9 n2/3 

^ ' 3 ^ 3' imu + msy^32E^k' ' 

<4o>^-k^)^^(:^)^/^ (13) 
3 3 TUs + nid 32E^k 

where the relative energy is calculated as E^. = (A/^ — to^ — m2)/2{'mi +1712) as mentioned previously from the energy 
E[ — (Af^ — 771^ — 777,2 )/2m2 rcst frame of 7772, and E[ is calculated from the Lorentz invariance of the momentum 

squared (pi + P2)c.m. — (Pi + P^^cst frame of ^^'^ energy calculated below in the center of mass frame. The 
index 1 of the two particles is assigned to the lighter one through out this letter. Since the number of the unknown 
parameters is more than that of the equations, more physical processes are needed to calculate. So the decay rates 
for the above meson are calculated in the following. 

From the relation dN = pvda, where N is the number of transition particles and v is the relative velocity between 
the two initial particles, the differential decay rate can be defined by 

S\m^....2,^ .4.4. ^ , , ^ rf^fci 'Pk2 

- 2^12^1-^1 (2") ^ +P2 - fci - fe^) (2^)32^, i2.f2u2 ' (^4) 

where S" is a symmetric factor for the initial particles and |5'(0)p = 'iErk/3'K from the Eq. ([III). The invariant 
amplitude for the pion decay process d{p2) + u{pi) — > p~(k2) + v{ki) is calculated as 

= 32G],cos^eH^,^U"', (15) 

where the hadronic tensor means that H^i^ = ^Tr[u[pi)^^^[l — 75)^(^2)^(^2)71^^(1 — 75)'^(pi) and the leptonic tensor 
is similar to that for the momentum — fci and fc2. The integration of the product of the tensors over the whole solid 
angle is calculated as 



j H^,L>"'dn = y 4pi • kip2 ■ k2dn 



„2^2 

^6A7TEiE2iujiLU2 + iy-Er), (16) 

0£/i£/2 

where the energies and momenta are calculated from the relation = A'P = {pi +^2)^ = (^i + ^2)^ in the center of 
mass frame as 

2lf^^^'+"''""'^), 



4M^ 
1 

2M~ 



^2^^iM^+ml), (17) 
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where M is the meson mass which is the bound system of mi and m2. From the above knowledge, the total decay 
rate is calculated as 



Stt ' ' ^ ^ 37rm2 3(1 _ _ ^)2') 

where the terms in the square bracket correspond to the pion decay constant in the usual calculation for 7r~ — > ^~+u 
and the kaon decay constant for K~ —> fj,~ + D. From the comparison between the decay and the charge radius 
of the pion in Refs. |^,|^, the space wave function |^'(0)p in the radius is 19 times larger than that in the decay. The 
reason is conjectured that the symmetric factor is not properly considered and the form factor of the pion is a dipole 
type in the reference. As the numerical calculations are shown below, the space wave function in the charge radii 
agrees to that in the decay rate strikingly, if the symmetric factor S is assumed to be 1/(4 x 9). The 1/4 means that 
the pseudoscalar meson is in the spin signet state, which may have been considered in Rcfs. [^[jsj already. However 
the color factor 1/9, which means that the decay takes place in the color singlet state and does not occur during the 
interactions of the exchange of 8 gluons, must have been ignored. 

For the typical current quark masses to„ = 10 MeV, = 28.61 MeV and nig = 160 MeV, the numerical values 
are calculated and compared with the experimental values px|-p^ in the second column of Table |, using the value 
k = 34655.5 MeV^ from the relation just above Eq. ( pi] ) with the estimated value as = 13.35 in Ref. |jl[. The reason 
for using the same value k through out all the physical quantities in the table is that the mass of the strange quark is 
known to be so heavy that it can not satisfy the condition 4m^ < q^. Instead of the decay rate, the decay constants 
are regarded as the experimental values for simplicity. The variation of the factor in the brace bracket in Eq. (|l^ ) 
is much small as the variation of the quark mass, actually it is less than 0.1 around its values 1.70 for the pion 
decay and 1.30 for the kaon decay. Hence the u- and d-quark mass ratio is extracted as 0.35 from the comparison of 
Eqs. (^3|) and (|l^ ) for the common E^k in the case of the pion, but the u- and s-quark mass ratio does not give an 
acceptable value in case of the kaon because of the small value of the kaon decay constant. If the calculated value 
of the kaon decay constant is replaced to the experimental value, the the u- and s-quark mass ratio is calculated 
as 0.065 and the CKM matrix element pKisI is estimated as 0.101. This matrix element is out of the range of the 
present standard values 0.217-0.224 ||l3| , [l^ . However there may be sufficient room to account for this discrepancy, 
if the Ke3 decay — > tt" -I- e"*" -I- is reconsidered carefully with an elegant quark model and the CKM matrix is 
chosen as the Kobayashi and Maskawa's original parameterization. This numerical calculation gives a new relation 

This comparison of the equations can be applied to other models, such as the wavefunction \E'(r) = ■^/a'^/Tre^'"' for 
the Coulomb potential Vcs — ~a/r and ^'(r) ~ i^aj-Kf'l^e^"''^ for the harmonic oscillator potential K;ff — br^. If 
the above quark masses are adjusted to the charge radii and taken as the inputs, the decay constants of the pion and 
kaon are 1.55 and 2.77 MeV for the coulomb potential and 274.2 and 583.9 MeV for the harmonic oscillator potential, 
respectively. Since the exponential function of the wavefunction reflects the asymptotic behavior of the potential, only 
the linear potential gives consistent results for the two different physical quantities, that is, charge radii and decay 
rates. 

Finally other current quark masses such as m„ = 1.5 MeV rrid — 4.29 MeV, and = 14.77 MeV give also 
interesting results as shown in the third column of table | for other suitable value k = 5057.66 MeV^ which means 
as — 34.0 in the point of view of Ref. as used previously. These quark masses are adjusted to the new relation 
/tt/Zk = M-^/Mk. Since the quark mass ratios are reproduced similar to that of the inputs, it is more reliable to 
adjust the quark masses to the experimental values for the charge radii rather than the decay constants. If the quark 
masses are adjusted to the decay constants and taken as inputs, then the estimated values for the charge radii of 
kaons are too large to accept them with ignoring the experimental values. It is also interesting to fit the quark masses 
to the experimental values for different values k for the pion and kaons, respectively, but it is still impossible to fit 
the experimental values to the charge radii and the decay constant of the kaons simultaneously with the same current 
quark masses due to the small kaon decay constant. This is the principal reason that there are four parameter, that 
is, rriu, rrid, nis, and fc, and five physical quantities, namely, < >, < >, < r'^o >, /^r, and /k, but it is 
difficult to fix the parameters with those physical quantities. 
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physical quantities 


calculated values I 


calculated values II 


experimental values 


<<- > 


0.438 fm^ 


0.433 fm^ 


0.43 ± 0.016 fm^ 


< 4- > 


0.347 fm^ 


0.227 fm^ 


0.34 ± 0.05 fm^ 


<4o > 


- 0.145 fm^ 


- 0.083 fm^ 


- 0.054 ± 0.026 fm^ 




129.38 MeV 


131.15 MeV 


130.7 ± 0.1 ± 0.36 MeV 


fK 


347.52 MeV 


463.89 MeV 


159.8 ± 1.4 ± 0.44 MeV 



TABLE I. The charge radius and the decay constant of the pion agree excellently to the experimental values, but those for 
the kaons do not account for the experimental values properly. The quark masses are adjusted to the radius of in the 
first calculation, and to the relation j-nj jK = Mt^/Mk, that is, the calculated /k in the second calculation. The experimental 
values are taken from Refs. [10-14]. 
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